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摘 要 : 针对 分 数 阶 传 里 叶 变 换 〈fractional fourier transform, FRFT) 对 Chirp 信号 进行 多 径 时 延 估 计时 的 不 足 ， 改 
进 了 一 种 按照 多 径 分 量 能 量 大 小 依次 消除 的 FRFT 多 径 时 延 估 计算 法 。 该 算法 以 选 代 方 式 进行 ， 按 多 径 分 量 的 能 量 
大 小 依次 返回 多 径 分 量 的 估计 值 ; 在 每 次 迭代 中 ， 包 含 子 迭代 以 准确 判定 当前 分 量 的 多 径 参 数 和 起 止 时 间 ， 然 后 利 
用 当前 多 径 参 数 生成 探测 信号 时 域 副 本 ， 将 其 从 残余 信号 中 减 去 ， 达 到 多 径 信号 的 分 离 。 以 根据 功率 时 延 分 布 的 特 
点 拟定 的 信道 模型 作为 传输 环境 ， 对 该 算法 进行 了 仿真 验证 。 仿 真 表明 ， 相 比 于 其 他 三 种 时 延 估 计算 法 ， 改 进 算 法 
能 够 更 准确 的 对 多 径 时 延 进 行 估计 。 
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Research on multipath delay estimation algorithm for chirp signal 


Feng Yali, Liu Cheng, Zhang Xiangbin 
(School of Communication & Information Engineering, Chongqing University of Posts & Telecommunications, Chongqing 
400065, China) 


Abstract: In order to solve the problem of the fractional Fourier transform (FRFT) in estimating the multipath delay for 
Chirp signal, this study proposed an improved FRFT multipath delay estimation algorithm based on the energy of each 
multipath component. The algorithm sequentially returns the estimated values of multipath components by iteration 
according to the energy magnitude of multipath components. In each iteration, there are a sub-iteration to accurately 
estimate the parameter and location of the current component, and then generates a time domain copy of the sounding signal, 
which is subtracted from the residual signal to reach separation of the multipath signal. This paper made some experiments 
which use channel model based on the characteristics of power delay distribution. and the results show that the improved 
algorithm can estimate the multipath delay more accurately than the other three delay estimation approaches. 

Key words: time delay estimation; chirp signal; fractional Fourier transform; multipath component cancellation 


0 引言 现 和 不 断 发 展 ， 作 为 非 平稳 信号 的 有 力 分 析 工 具 ， 因 其 变换 
di 基底 为 一 组 正 交 的 Chir 信号 ， 且 具有 多 种 与 快速 傅 里 叶 变 

Chirp 信号 又 称 线性 调频 信号 ， 存 在 于 自然 界 的 很 多 领 。” 换 (fast fourier tranform, FFT) 计算 相当 的 离散 算法 ， 近 年 来 
域 ， 如 天 体 物理 学 中 的 引力 波 、 蝙 蝠 的 回声 定位 、 医 学 中 的 “被 广泛 应 用 于 Chirp 信号 的 检测 与 参数 估计 人 中。Ozaktas 等 
音信 号 等 。 因 Chirp 信号 具有 高 处 理 增益 、 低 发 射 功率 、 人 中 于 1996 年 提出 了 一 种 精确 计算 FRFT 的 离散 算法 。 所 提 
抗 多 普 勒 频 移 等 优点 ， 所 以 其 在 雷达 帆 、 扩 频 通 信忠 、 水 下 出 的 算法 的 时 间 复 杂 度 为 O(NlogN )。 至 此 ，FRFT 不 再 停 
IUBIM UWB 探测 外 等 领域 中 有 着 相当 广泛 的 应 用 。 对  ” 留 于 理论 研究 阶段 ， 迈 向 了 工程 实践 的 道路 。Tao 等 人 中 在 
Chirp 信号 进行 多 径 时 延 估 计 是 信道 探测 中 对 测量 数据 进 
处 理 的 一 个 关键 环节 。 传 统 的 多 径 时 延 估计 主要 采用 匹配 滤 
、 解 卷 积 (de-convolution，DC) 算法 。Bell 等 人 提出 的 
滤波 算法 ， 其 时 延 分 辩 率 为 脉冲 压缩 后 sinc 函数 的 主办 
意 〈 近 似 于 探测 信号 带宽 的 倒数 )，sinc 函数 旁 兴 对 多 径 


S 
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前 人 的 研究 基础 上 ， 首 次 推导 了 FEFRFT 时 延 估 计 (FRFT time 
delay estimation,FRTDE ) 算 法 的 精度 与 SNR. 之 间 关 系 的 数学 
表达 式 。Zhou 等 人 09 基 于 FRITDE 算法 提出 了 一 种 新 的 子 样 
本 FRFT 时 间 延 迟 估计 (subsample FRFT time delay estimation, 
SFRTDE) 算 法 ， 它 进一步 利用 信号 FRFT 输出 峰值 的 相位 
息 的 提取 影响 比较 大 外 。 DC 算法 对 噪声 的 存在 较为 敏感 ， 息 来 搜索 匹配 的 时 延 ， 并 通过 计算 机 仿真 验证 了 其 相对 于 
有 在 信 噪 比 〈signal to noise ratio» SNR) 足够 好 的 时 候 才 FRTDE 和 平均 幅度 平方 分 数 差分 函数 算法 的 有 效 性 和 优越 
很 好 地 估计 出 多 径 时 延 名 ;传统 的 这 些 算 法 要 求 探测 信和 号 为 性 。Sejdié 等 人 HH 对 FRFT 做 了 一 个 较为 详细 的 概述 ， 有 具体 
平稳 信号 , M Chirp 信号 为 非 平稳 信号 ,所 以 用 其 来 处 理 Chirp 介绍 了 FRET 定义 ， 数 学 实现 及 其 应 用 领域 ， 将 其 和 数学 变 
言 号 存在 先天 不 足 。 随 后 的 一 些 基 于 高 阶 累 积 量 的 算法 ， 如 换 联系 起 来 ， 强 调 实现 FRFT 的 各 种 方法 ， 并 比较 了 它们 的 
互 四 阶 矩 最 小 范 数 时 延 估 计算 法 和 四 二 阶 归 一 化 累积 量 自 适 优 缺 点 ， 总 结 了 FRFT 处 理 非 平 稳 信 号 和 时 变 信 号 所 存在 的 
应 时 延 估 计 方 法 ， 此 类 算法 首先 把 Chirp 信号 转换 为 平稳 信 优势 。 
号 ， 再 进行 时 延 估 计 ， 时 延 估 计 准 确 度 高 ， 但 此 类 算法 时 间 在 信道 测量 系统 硬件 资源 和 计算 复杂 度 结合 考虑 的 情况 
复杂 度 普遍 较 大 ， 不 适合 实际 使 用 趾 。 随 着 FRFT 理论 的 出 下 ， 经 济 适用 的 多 径 时 延 估 计 方 法 成 为 当下 迫切 的 需求 。 本 
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文 从 现 有 的 时 延 估计 技术 入 手 , 同时 应 用 Chirp 信号 在 FRFT ”对 接收 到 的 带 有 了 噪声 的 多 径 信 号 r( 进行 分 数 阶 傅 里 叶 变换 ， 
域 的 峰值 汇聚 特性 和 时 域 的 脉冲 压缩 特性 来 研究 新 的 时 延 估 ”将 其 变换 到 使 其 能 量 产 生 聚 集 的 最 佳 阶 傅 里 叶 变 换 域 上 ， 然 
计算 法 ， 以 便 从 测量 数据 中 得 到 更 精确 的 信道 模型 。 后 在 分 数 域 上 对 信号 幅度 谱 的 峰值 位 置 进行 检测 ， 获 得 峰值 
1 ”信号 模型 和 峰值 点 位 置 的 坐标 ， 并 与 发 射 信 号 sO 在 分 数 域 上 幅度 谱 
的 峰值 位 置 进行 比较 ， 就 可 以 获得 幅度 谱 的 位 置 偏 移 量 及 峰 


Chirp 信道 探测 信号 的 时 域 表达 式 为 值 比 例 。 为 方便 推导 ,假设 M = 二 1， 幅 度 4 三 1， 根 据 分 数 阶 
SE -exp(2z fit jak?) -T,12<t<T,12 (1) — 傅 里 叶 变 换 的 时 移 特性 ， 则 有 
其 中 : 为 中 心 频率 ; 为 调频 率 ; TJ Chirp 信号 的 持续 F|s(t — 7)] 2 e/ t Psneemae ursing f (u—r cos) (6) 
时 间 ，Chirp 信号 的 带宽 B- KT, 。 探 测 信号 在 传播 过 程 中 可 能 对 两 边 取 绝 对 值得 
会 遇 到 建筑 物 、 树 木 ， 以 及 起 伏 的 地 形 ,这 就 会 导致 电波 出 现 IF-F,- reos D 
反射 、 散 射 、 透 射 及 绕 射 等 现象 ， 所 以 到 达 接 收 端 天 线 的 信 > 


号 不 是 只 从 单一 路 径 来 的 ,而 是 许多 路 径 的 信号 合成 。 多 径 信 
道 时 延 估计 问题 的 接收 信号 模型 可 以 表示 为 


r=) ast -T„)+ elt) 2) 


LH: e(t) 为 零 均 值 、 方差 为 o? 的 高 斯 白 噪 声 ; M x Ams f, 
为 未 知 参数 ， 分 别 对 应 信道 的 多 径 信 号 数目 ， 第 m 条 多 径 信 
号 的 幅度 和 时 延 。 

2 基本 的 FRFT 多 径 时 延 估 计算 法 

2.1 Chirp 信号 的 FRFT 


归 一 化 幅度 


100 


x. pi — Vus TD cu 变换 阶 次 0 分 数 域 采样 点 
Chirp 信号 是 FRFT 的 标准 正 交 基 。 根 据 正 交 基 的 特性 ， : dicis 
Chirp 信号 变换 到 某 个 特定 阶 次 的 分 数 阶 傅 里 时 变换 域 上 将 图 1 各 阶 次 的 FRFT 波形 
是 一 个 冲击 函数 ， 产 生 能 量 聚 集 的 特性 ， 所 以 FRFT 特别 适 Fig.1 FRFT waveform of each order 
合 处 理 Chirp 信号 中, 实际 应 用 中 Chirp 信号 都 是 有 限 长 的 ， 1 
即 Chirp 信号 的 持续 时 间 是 有 限 的 ， 因 此 其 分 数 阶 传 里 叶 变 os- 
换 的 形式 也 有 所 改变 。Chirp 探测 信号 P Bri ERRA os. 
F"[s(t)] 7 f, (0) 2 J1— jcota x 07r 
a dp 0.6 F 
| exp[ jx(t? cota — 2ut csc æ +u? cota)]* E 
-T2 3) p^ 
exp[j(2 f.t + zkt?)]dt = 41 — jcoto -exp( jru? cota) x ( E our 
i exp[jzt? (cotæ - k)]- exp[ 2ztt( f, —u csc axo) ]dt dl 
-n2 0.2 
其 中 : a 为 旋转 角度 ，a= pz/2 。 最 佳 变 换 阶 次 与 Chirp 信号 zu d 
Vis AIR Z8 Je Us o 20 40 60 80 100 120 140 160 180 200 
1 分 数 域 采 样 点 
PRA a 图 2 最 佳 阶 次 的 FRFT 波形 
zn Fig.2  FRFT waveform of optimal order 
即 当 k=-cota 时 , Chirp 信号 变换 到 对 应 的 分 数 阶 传 里 叶 变 换 HR (7) 可 以 看 出 , 信号 在 时 域 上 的 时 延 变 换 到 分 数 域 
域 表 现 出 sinc 函数 的 形式 ， 式 (4) 可 化 简 为 上 表现 为 信号 幅度 谱 平 移 rcosa 。 只 要 获得 时 延 信号 在 分 数 
1,0) = JI Joora exp jt cota): 域 上 幅度 谱 的 位 置 偏 移 量 就 可 以 获得 对 时 延 的 估计 ， 
Ls t - |i, -û,|/ cosa (8) 
| exp[ j2zt( f. -ucsca)]dt = (5) 
Edd trh, Â -agmax(F^GG) ， 为 探测 信号 的 峰值 位 
T, 41— jcota :exp(jAxu? cota):sinc[xT.(f. —ucsca)] s " 为 探测 信 号 的 i Li 位 
?4 f, -ucsca , sinc 函数 取 最 大 值 为 1， 则 此 时 1/5 G0| Bo i= P : -— i " 
~o d cemme^00) ， 为 接收 信号 的 峰值 位 置 。 根 据 分 数 阶 伟 


最 大 值 . 从 以 上 推导 过 程 可 以 看 出 , 当 FRET 的 旋转 角度 o 满 


足 时 上 =-cota 的 条 件 时 , Chirp 信号 在 频 域 出 现 了 能 量 高 度 江 里 叶 变换 的 线性 特性 ， 则 

聚 的 特性 ， 被 压缩 成 了 一 条 “ 谱 线 ”， 包 含 了 Chirp 信号 的 全 F"[as(t —7)] s aF^[s(t - 7)] - 0) 
部 能 量 ， 即 Chirp 的 频 域 压缩 。Chirp 信号 各 阶 次 的 分 数 域 三 ae) Posee sne f, (u— reos an) 

维 波形 及 最 佳 阶 次 的 分 数 域 波形 分 别 如 图 1 和 2 所 示 。 对 两 边 取 绝对 值 ， 得 


2.2 FRTDE 算法 

FRTDE 算法 主要 利用 Chirp 信号 在 FRFT 域 的 峰值 汇聚 
特性 对 多 径 时 延 进 行 估计 。Chirp 信号 在 分 数 阶 傅 里 叶 变 换 同村 
域 上 有 良好 的 能 量 聚 集 特 性 ， 而 高 斯 白 噪声 变换 到 分 数 域 仍 数 域 上 可 
为 高 斯 白 噪 声 ， 不 会 在 任何 阶 分 数 域 上 产生 能 量 聚 集 ， 可 以 RELE 
利用 Chirp 信号 和 高 斯 白 噪 声 的 此 特性 进行 多 径 时 延 估计 [09。 的 比例 ， 


|F"[as(t — 7)]| laf, (u cosa) alf, (u zcosaw)| (10) 


A (10) 可 以 看 出 ， 多 径 信号 的 幅 值 变换 到 分 
于 在 探测 信号 的 分 数 域 峰值 上 乘 以 其 幅 值 ， 只 要 
号 在 分 数 域 上 的 峰值 和 探测 信号 在 分 数 域 的 峰值 
可 以 获得 对 多 径 信号 幅度 的 估计 : 


S up E 
di 
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~ max(F^(r)) 
“= max(F? GG) 00 


假设 接收 到 的 信号 为 来 自 单一 路 径 的 信号 ， 利 用 基本 的 
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MCEE-FRTDE 以 迭代 方式 进行 , 按 多 径 分 量 的 能 量 大 小 依次 
返回 多 径 分 量 的 估计 值 ， 每 次 返回 多 径 分 量 估计 值 时 ， 对 该 
分 量 进 行 重建 ， 同 时 在 残余 信号 减 去 当前 分 量 以 更 新 残余 信 


分 数 阶 傅 里 叶 变 换算 法 可 以 较为 精确 地 对 多 径 时 延 和 幅度 进 
行 估计 。 然 而 对 于 实际 的 传播 环境 ， 接 收 信号 为 来 自 各 条 路 
径 的 多 分 量 信 号 的 合 加 ， 各 路 径 分 量 会 不 可 避免 的 对 彼此 产 
生 影 响 ， 基 本 的 分 数 阶 傅 里 叶 变 换 时 延 估 计算 法 未 考虑 到 路 
径 之 间 的 相互 影响 ,信道 估计 的 精度 会 受到 路 径 分 量 的 影响 。 


3 ”改进 的 FRFT 多 径 时 延 估 计算 法 


经 过 多 径 信道 的 传播 后 ， 接 收 信号 为 来 自 各 条 路 径 的 多 
分 量 信号 的 个 加 ,各 路 径 分 量 会 不 可 避免 的 对 彼此 产生 影响 。 


mi 


; WERE SUE TY 982r fei m IU. ARERR 
牛 为 残余 信号 的 最 佳 阶 次 的 FRFT 不 再 具有 尖峰 特性 ， 此 理 
论 依据 为 2.2 节 中 提 到 的 Chirp 信号 和 高 斯 白 噪声 在 FRFT 域 
的 特性 ， 同 理 ， 脉 冲 干扰 在 频谱 的 特征 类 似 于 高 斯 白 噪声 ， 
也 不 会 在 任何 阶 次 产生 能 量 汇聚 。 在 每 次 迭代 中 ， 包 含 子 迭 
代 以 准确 的 判定 当前 分 量 的 多 径 参 数 和 其 在 残余 信号 中 的 具 
体位 置 。 多 径 分 量 和 残余 信号 中 的 示意 图 如 图 3 所 示 。 在 子 
RF, ME Chirp 信号 的 匹配 滤波 特性 来 找到 当前 分 量 在 
余 信 号 中 所 处 的 具体 位 置 ， 然 后 对 当前 分 量 所 处 的 信号 段 


a 


受信 号 分 离 理论 思想 的 启发 ， 在 信道 估计 中 ， 可 以 对 各 路 径 
信和 号 按照 能 量 大 小 进行 依次 剥离 进而 消除 高 能 量 信号 对 其 他 
路 径 分 量 估计 的 影响 ， 本 文 改进 了 一 种 按照 多 径 分 量 能 量 大 
小 依次 消除 的 FRTDE Cmultipath components are eliminated in 
order of energy-based FRTDE,MCEE-FRTDE). 

接收 信号 rO E M E E [Tsar Tena] 上 包含 第 mn 条 多 径 的 
持续 时 间 [Tmstar,Tmend]。 假设 接收 信号 中 第 条 信号 分 量 的 幅 
度 w 最 大 ， 式 OO 可 改写 为 


rOD=wO+ > set) (12) 


3.1 多 径 分 量 消除 方法 

多 径 分 量 消除 方法 的 思想 是 : 首先 通过 FRFT 算法 估计 
出 接收 信号 中 能 量 最 强 的 信号 的 时 延 和 幅度 ， 重 建 该 多 径 分 
量 ， 然 后 根据 Chirp 信号 的 匹配 滤波 特性 找 出 该 分 量 在 接收 
信号 中 的 具体 位 置 ， 最 后 从 接收 信号 的 对 应 位 置 将 其 减 去 ， 
直至 接收 信号 中 无 多 径 分 量 为 止 。 
基于 上 一 节 的 讨论 , 可 以 通过 式 (13) 和 (14) 求 出 第 
条 分 量 信 号 对 应 的 时 延 和 幅度 : 


liz 
jury 


£, «|i, - à,|/ cos æ (13) 
,  max(F*(r()) 
” max(F^ (s(t) (14) 
根据 f, 和 à, 重建 第 n 条 信 号 分 量 s,Q) : 
s, (0) 2 à,s(t —2,) (15) 
根据 Chirp 信号 的 匹配 滤波 特性 ，w (0 KRIE s rO 


的 匹配 输出 为 一 个 脉冲 压缩 信号 ， 由 脉冲 压缩 信号 的 最 大 值 
所 在 位 置 可 求 出 sO Æ rO 中 存在 的 具体 位 置 。 


万 (1) =r(t)*s*,(t)= i Sn(t—T)s*,(T)dT+ 
"um 7 M (16) 
>» f 5-05, (Ddr+ | e(t—7)s", (T)dr 


类 似 于 文献 [15], 上 式 右 边 第 一 项 为 匹配 输出 ， 可 表示 为 


oo 


E (t—7)s;(r)dr = a 7) 


T A((t — 7, )/T)sin c[zKT (t — 7, )A((t — T, )/T)]e 7e) 

中 : A(O 为 定义 在 [-1,1] 上 的 正三 角 脉 冲 函数 , 在 [-1,0] 上 从 0 
线性 增长 到 1， 在 [0,1] 上 从 1 线性 减 小 到 0。 式 (16) 中 的 第 
二 项 为 能 量 较 弱 的 分 量 的 非 匹 配 输出 ， 视 为 干扰 。 

结合 Chirp 信号 在 FRFT 域 的 峰值 汇聚 特性 和 时 域 的 脉 
冲压 缩 特性 , 能 够 从 式 (12) 中 将 分 量 s OO 与 其 他 分 量 分 离 ， 
在 式 (16) 中 的 第 一 项 即 被 移 除 , 在 分 离 能 量 次 强 的 分 量 时 ， 
ERAT w(0D 残余 能 量 的 干扰 。 
3.2 MCEE-FRTDE 算法 

MCEE-FRTDE 算法 将 同时 应 用 Chirp 信号 在 FRFT 域 的 
峰值 汇聚 特性 和 时 域 的 脉冲 压缩 特性 对 多 径 时 延 进行 估计 。 


AL 


再 次 做 FRFT 以 估计 当前 分 量 参数 ， 友 代 结 束 条 件 为 两 次 估 
计 的 时 延 和 幅度 相差 小 于 一 个 极 小 量 ， 即 认为 估计 值 稳定 。 

残余 信号 总 时 长 
[Istart.Tend] 


t 


Tnstart Tnend Tend 


当前 分 量 时 长 
[Tnstart .Tnendl 
图 3 当前 分 量 与 残余 信号 示意 图 
Fig.3 Diagram of current component and residual signal 
基于 以 上 的 分 析 ，MCEE-FRTDE 算法 具体 步骤 如 下 ; 
输入 : 探测 信号 sO ; 接收 信号 rO ; 可 接受 误差 小 量 a， 


£y 

初始 化 : 残余 信号 nO Tr. RRI n=1; 

a) 分 别 对 DO RU nO) 做 最 佳 阶 传 里 叶 变 换 
(p= Sy 


b) 通过 式 (13) 和 (14) 分 别 求 出 第 个 分 量 信号 的 名 
fll à, 。 

c) HE S. a WX (15) 重建 信号 w(D) 。 

d) 计算 sO Bd gu [9 5 nO 的 卷 积 、s,(0 Bodegas s 
与 自身 的 卷 积 ,得 到 w() 在 i() 中 存在 的 具体 位 置 , 记 为 i。 


i-es 7 as] (18) 

其 中 
isa, = arg max (s; (t) * s, (1)) (19) 
im = arg max(s; (t) * r, (0) (20) 


e) 重新 对 搜索 出 来 的 信号 段 ni T1 f dec P Y XCTI 
FRFT， 获 得 一 组 新 的 名 和 o 

D 令 ar, = 各 一 名 | , aa, =â, -à,| o 如果 Ar > 6 或 者 sa >&， 
4 £26, à-à,, JP IbURGC)-e) 否则 ， 进 入 下 一 步 。 

g) 更 新 残余 信号 的 信号 段 i[i:i+7]: 

nli:i+T]=r[i:i+T]-s,(t) (21) 

h) o nO 做 最 佳 阶 次 的 FRFT， 如 果 不 再 有 尖峰 ， 即 多 
径 分 量 分 离 完毕 ， 则 停止 迭代 ; 否则 ， 令 n=n+1， 从 步骤 b) 
开始 重复 以 上 步骤 。 


输出 : 多 径 时 延 参数 估计 值 Ue 


从 MCEE-FRTDE 算法 的 具体 步骤 可 知 ， 其 时 延 分 辩 率 
主要 取决 于 Chirp 信号 匹配 滤波 的 时 延 分 辩 率 和 最 佳 阶 傅 里 
叶 变 换 的 时 延 分 辩 率 ,对 多 径 分 量 进行 消除 的 目的 是 减少 干 
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扰 ， 使 得 时 延 估计 结果 更 加 准确 。 由 文献 [16] 可 知 ， 匹 配 滤 杂 度 ， 一 次 卷 积 计算 需要 两 次 FFT、 一 次 IFFT 和 一 次 乘法 


波 和 最 佳 阶 传 里 叶 变 换 的 时 延 分 辩 率 都 约 为 WB， 所 以 ” 运算 , 所 以 步骤 g 四 的 复杂 度 为 230(NlogN)+0O(N))=0O(NlogN); 

MEBESERTDE SA 步骤 6) 和 外 需要 进行 一 次 FRFT， 复 杂 度 O(NlogN) 。 综 上 ， 
2i MCEE-FRTDE 算法 的 时 i 杂 度 为 OOVlogN) 。 DC 算法 需 

4 ”仿真 结果 及 性 能 分 析 算法 的 时 间 复杂 度 为 O(N1og 算法 需要 


两 次 FFT、 一 次 IFFT 和 一 次 除法 运算 , 然后 进行 峰值 搜索 ; 
Chirp 探测 信号 的 仿真 参数 如 表 1 所 示 。 功 率 时 延 分 布 ”FRTDE 和 SFRTDE 均 只 需要 进行 两 次 FRFT， 然后 进行 峰值 
一 般 服 从 指数 分 布 01， 可 通过 这 一 性 质 拟 定 信道 模型 参数 ， 搜索 , 所 以 解 卷 积 算法 、FRTDE 和 SFRTDE 的 算法 复杂 度 也 
根据 功率 时 延 分 布 的 特点 拟定 的 信道 模型 参数 ( 表 2), 过 采 都 为 COVlogN) ， 与 MCEE-FRTDE 算法 的 复杂 度 是 相当 的 。 
样 因 子 为 5。 每 次 仿真 都 在 信道 上 附加 一 定 SNR 的 高 斯 白 噪 “虽然 上 述 各 种 算法 的 时 间 复 杂 度 相当 ， 但 是 因为 
声 和 随机 脉冲 干扰 ， 脉 冲 干扰 出 现 的 概率 为 0.1， 信 噪 比 服 ^^ MCEE-FRTDE 存在 迭代 以 优化 估计 出 来 的 参数 值 , 其 仿真 运 
从 [-20dB, 0dB] 之 间 的 均匀 分 布 。 采 用 500 次 Monte Carlo fj; ” 行 时 间 高 于 其 他 算法 ， 所 以 其 在 运行 效率 方面 仍 有 改进 的 空 
真 的 平均 来 得 到 多 径 信 道 时 延 估计 的 归 一 化 均 方 误差 ” 间 。 
(normalized mean square error, NMSE)。NMSE 由 式 (22) 


N 


计算 。 | | wem 
0.08 F —4— SFRTDE 
M 2 一 此 一 FRTDE 
937 一 名 | 0.07 上 全 De 
NMSE ==] (22) 
A : 
表 1 Chip 探测 信号 的 参数 ii 
Table 1 Parameters of Chirp detection signal 2 
参数 值 
T; 1.024 us 
fe 490 MHz | == 
-15 -10 -5 0 5 
B 100 MHz SNR(dB) 
表 2 拟定 的 信道 模型 参数 图 4 ”时 延 估 计 的 归 一 化 均 方 误差 与 信 品 比 的 关系 
Table 2 Parameters of proposed channel model Fig.4 Relationship between NMSE of time delay estimation and 
路 径 时 延 /ns 幅度 /dB SNR 
1 0 0 2 l | 
- -1 s Mil Wl 
3 100 33 UN WW. ANM 
4 170 -5.6 ^o o 26 RNC so ex 70 
5 240 -7.9 5 T r r r EE 
为 了 更 好 地 对 改进 算法 的 性 能 作出 评价 ， 本 文 将 E «Mode Phot An 
MCEE-FRTDE 算法 和 DC 算法 、FRTDE 算法 以 及 SFRTDE 
算法 进行 分 析 比 较 。 图 4 为 几 种 算法 估计 时 延 的 归 一 化 均 方 BEER s 
误差 与 信 噪 比 的 关系 。 总 体 而 言 ， 各 种 算法 的 性 能 都 随 信 噪 NINE qum 
比 的 增加 而 变 得 越 来 越 好 。 在 图 4 中 , DC 算法 的 性 能 最 差 ， DW | 
A DC 算法 对 噪声 较为 敏感 ， 特 别 是 在 低 信 噪 比 的 时 候 ， 7^» we — me oo A9 SE sb Too 
噪声 强度 较 大 ， 其 估计 准确 度 受 到 较 大 影响 ，SFRIDE 算法 xs 
利用 信号 FRFT 输出 峰值 的 相位 信息 来 搜索 匹配 的 时 延 ， 比 J5 ”多 尺度 多 时 延 源 信号 的 恢复 
FRTDE 算法 直接 利用 峰值 信息 具有 更 高 的 时 延 估计 精度 ; Fig.5 Recovery of multi-scale and multi-delay source signal 
MCEE-FRTDE 算法 考虑 到 了 各 路 径 分 量 之 间 的 相互 影响 , 并 


5 ”结束 语 


日 还 通过 迭代 使 得 FRFT 长 度 能 够 自 适应 地 调整 ， 最 终 达 到 
目标 分 量 的 起 始 位 置 。 相 比 以 上 算法 , MCEE-FRTDE 算法 能 本 文 改 进 了 一 种 按照 多 径 分 量 能 量 大 小 依次 消除 的 新 方 
够 更 准确 对 时 延 进行 估计 。 法 MCEE-FRTDE。 在 该 方法 中 ， 利 用 FRFT 对 Chirp 信号 的 


利用 MCEE-FRTDE 算法 估计 得 到 的 多 径 信 道 参 数 ， 可 ”汇聚 特性 ， 以 迭代 的 方式 ， 按 照 路 径 能 量 从 大 到 小 的 顺序 将 
以 将 多 尺度 多 时 延 的 源 信号 恢复 出 来 。 如 图 5 Bras. 上 到 ”多 径 分 量 的 位 置 及 参数 估计 出 来 ， 并 从 接收 信号 中 分 离 。 与 
下 显示 的 分 别 是 源 信号 、 加 噪声 (SNR=1dB) 的 接收 信号 ， DC, FRTDE 和 SFRTDE 相 比 , 所 提出 的 时 延 估计 方法 可 以 为 
以 及 用 多 径 信道 参数 估计 值 重 构 出 来 的 源 信号 与 源 信号 的 比 Chip 信号 提供 更 好 的 时 延 估计 结果 ， 但 同时 也 增加 了 测量 


较 。 可 以 看 到 重 构 的 源 信号 与 源 信号 比较 吻合 ， 再 次 说 明了 ”系统 对 其 进行 处 理 的 时 间 。 在 实际 进行 信道 测量 的 过 程 中 ， 
改进 算法 的 准确 性 。 一 般 都 是 对 接收 端 信息 进行 离线 处 理 ， 对 算法 的 实时 性 要 求 

为 了 评估 MCEE-FRTDE 算法 的 计算 复杂 度 和 实用 性 ， 不 高 。 综 合算 法 的 估计 精度 、 时 间 复 杂 度 和 系统 要 求 来 说 ， 
需要 对 算法 进行 时 间 复 杂 度 的 对 比 。 假 设 信号 向 量 长 度 都 为 MCEE-FRTDE 算法 具有 很 好 的 性 能 。 


杂 度 为 20(NlogN)=O(NlogN) ; 步骤 由 需要 进行 两 次 卷 积 运算 ， I 
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